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We have investigated the isoscalar giant monopole resonances (ISGMR) 
in 112 ~ 124 Sn and 106 ~ 116 Cd nuclei using inelastic scattering of 386-MeV 
a-particles at extremely forward angles, including 0°. The strength distri- 
butions for various multipoles were extracted by a multipole decomposition 
analysis based on the expected angular distributions of the respective mul- 
tipoles. From the ISGMR results, a value of K T ~ —500 MeV is obtained 
for the asymmetry term in the nuclear incomprcssibility. 

PACS numbers: 24.30.Cz; 21.65.+f; 25.55.Ci; 27.40. +z 

1. Introduction 

Investigation of giant resonances has had a rich history over the past six 
decades since the discovery of the isovector giant dipole resonance (IVGDR) [lj 
[21 [3] . Of the various giant resonances that have been identified and inves- 
tigated since, the compression-mode giant resonances-the isoscalar giant 
monopole resonance (ISGMR; the "breathing mode") and the isoscalar gi- 
ant dipole resonance (ISGDR; the "squeezing mode" )-occupy the pride of 
place because their energy is directly related to the nuclear incompressibil- 
ity, Ka, from which the incompressibility of infinite nuclear matter, K^, 
may be deduced [I] . The latter is critical in our understanding of a number 
of interesting phenomena from collective excitations of nuclei to supernova 
explosions and radii of neutron stars. 

In this review, we discuss some recent results on the ISGMR over a 
series of isotopes of Sn and Cd. This data has provided an "experimental" 
value for the asymmetry term, K T , of nuclear incompressibility; this term, 
associated with the neutron excess (N — Z) is important, for example, in 
obtaining radii of neutron stars in EOS calculations. 
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2. Experimental Techniques 

The experiment was performed at the ring cyclotron facility of the Re- 
search Center for Nuclear Physics (RCNP), Osaka University, using inelas- 
tic scattering of 386- MeV a particles at extremely forward angles, includ- 
ing 0°. Details of the experimental procedures and data- analysis tech- 
niques have been provided elsewhere [6J and are only briefly described 
here. Self-supporting target foils of enriched 112 . 114 > n 6,ii8, 120,122, i24g n anc j 
106,110,i12,114 ) i16q ( j i so t pes of thickness ~5-10 mg/cm 2 were employed; spe- 
cial target frames with a large aperture were used in order to reduce the 
background caused by the beam-halo hitting the frames. Data were also 
taken with a nat C target at the actual field settings used in the experiments 
and energy calibration was obtained from the peak positions of the 7.652- 
and 9.641-MeV states in the 12 C(a,a') spectra. 

Inelastically-scattered a particles were momentum-analyzed with the 
high-resolution magnetic spectrometer "Grand Raiden" [7] and the vertical 
and horizontal positions of the a particles were measured with a focal-plane 
detector system comprised of two position-sensitive multi-wire drift cham- 
bers (MWDCs) and two scintillators [8]. The MWDCs and scintillators en- 
abled us to make particle identification and to reconstruct the trajectories of 
the scattered particles. The scattering angle at the target and the momen- 
tum of the scattered particles were determined by the ray-tracing method. 
The vertical-position spectrum obtained in the double-focusing mode of the 
spectrometer was exploited to eliminate the instrumental background [8, 9j. 
Examples of such "free-from-instrumental-background" inelastic scattering 
spectra for the Sn isotopes have been provided in Refs. [6]. 

A multipole decomposition analysis (MDA) procedure [10] was employed 
to extract the strengths of the ISGMR, along with those for other isoscalar 
resonances (up to L = 3). The associated DWBA calculations were per- 
formed with the computer code PTOLEMY [11] . following the method of 
Satchler and Khoa [12] . using the density-dependent single- folding model 
for the real part, obtained with a Gaussian a-nucleon potential, and a phe- 
nomenological Woods-Saxon potential for the imaginary term. The optical- 
model (OM) parameters were determined by fitting the differential cross 
sections of elastic a scattering measured in a companion experiment; the 
efficacy of the OM parameters was tested by comparing the experimental 
cross sections for the first 2 + states in these nuclei with calculated cross sec- 
tions using a collective form factor and previously-established B{E2) values. 
Again, examples of MDA fits to the experimental angular distributions of 
the differential cross sections have been provided in Refs. [BJ. 
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3. Results and Discussion 

We have extracted strength distributions for L=0, 1, 2, and 3 multipoles 
over the energy range 8.5 MeV-31.5 MeV in all the Sn and Cd isotopes 
investigated in this work. The ISGMR strength distributions for the Cd 
isotopes are presented in Fig. 1. The L = strength distributions were 
fitted with a Lorentzian function to determine the centroid energies and 
widths of the ISGMR. These fits are shown superimposed in Fig. [TJ the 
corresponding fitting parameters are presented in Table 1. Also presented 
are the various moment ratios for the experimental ISGMR strength distri- 
butions calculated over the excitation-energy range, E x = 10.5-20.5 MeV, 
encompassing the ISGMR peak. The results for the Sn isotopes have been 
presented previously [5l [6] . 
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Fig. 1. ISGMR strength distributions obtained for the Cd isotopes in the present 
experiment. Error bars represent the uncertainties from fitting the angular distri- 
butions in the MDA procedure. The solid lines show Lorentzian fits to the data. 
These results should be considered preliminary. 

A note of caution is apt here. There is a small, near-constant ISGMR 
strength up to the highest excitation energies measured in this experiment. 
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The exact reason behind this extra strength is not quite well understood. 
However, similarly enhanced El strengths at high excitation energies were 
noted previously [13] and have been attributed to contributions to the 
continuum from three-body channels, such as knockout reactions [13] • These 
processes are implicitly included in the MDA as background and may lead 
to spurious contributions to the extracted multipole strengths at higher 
energies where the associated cross sections are very small. 

Table 1. Lorentzian-fit parameters, as extracted from MDA and the various mo- 
ment ratios for the ISGMR strength distributions in the Cd isotopes. All moments 
have been calculated over E x = 10.5-20.5 MeV. These results should be con- 
sidered preliminary. 



Target E GMR (MeV) F(MeV) g (MeV) (MeV) (MeV) 



106 Cd 


16.4 ±0.1 


4.4 ±0.5 


16.7 ±0.1 


16.5 ±0.2 


17.2 ±0.3 


110 Cd 


16.1 ±0.1 


4.1 ±0.6 


16.3 ±0.1 


16.1 ±0.2 


16.9 ±0.3 


112 Cd 


15.8 ±0.1 


4.9 ±0.7 


16.2 ±0.1 


16.0 ±0.2 


16.8 ±0.2 


114 Cd 


15.5 ±0.1 


5.0 ±0.6 


16.1 ±0.1 


15.8 ±0.2 


16.7 ±0.4 


116 Cd 


15.4 ±0.1 


5.0 ±0.4 


15.9 ±0.1 


15.7 ±0.1 


16.6 ±0.3 



The excitation energy of the ISGMR is expressed in the scaling model |15] 

as: 



where m is the nucleon mass, < r 2 > the ground-state mean-square radius, 
and Ka, the incompressibility of the nucleus. 

Further, the incompressibility of a nucleus, Ka, may be expressed as: 



K A ~ K vol {\ + cA- 1 / 3 ) + K T {{N - Z)/A) 2 + K Coul Z 2 A~ A / 3 (2) 



Here, c « -1 (see, for example, Ref. [16], and KqouI is essentially model 
independent (in the sense that the deviations from one theoretical model 
to another are quite small), so that the associated term can be calcu- 
lated for a given isotope. Thus, for a series of isotopes, the difference 
Ka — KcouiZ 2 A~ i / 3 may be approximated to have a quadratic relation- 
ship with the asymmetry parameter ((N - Z)/A)), of the type y = A + 
Bx 2 , with K T being the coefficient, B, of the quadratic term. It has been 
established previously [T71 [18] that direct fits to the Eq. 2 do not provide 
good constraints on the value of K^. However, this expression has been 
used in this case not to obtain a value for K^, but only to demonstrate the 
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approximately quadratic relationship between Ka and the asymmetry pa- 
rameter. In addition, it should be understood that this expression provides 
only an "average" value for K T and the available data is not sensitive to 
higher-order effects like the "surface" part of this term. Here, c ~ -1 (see, 
for example, Ref. [15] . and Kc ou i is essentially model independent (in the 
sense that the deviations from one theoretical model to another are quite 
small), so that the associated term can be calculated for a given isotope. 
Thus, for a series of isotopes, the difference Ka — Kc uiZ 2 A~ A ^ may be 
approximated to have a quadratic relationship with the asymmetry param- 
eter ((N - Z)/A)), of the type y = A + Bx 2 , with K T being the coefficient, 
B, of the quadratic term. It has been established previously [TTJ, [T2] that 
direct fits to the Eq. 2 do not provide good constraints on the value of K^. 
However, this expression has been used in this case not to obtain a value 
for .ftTooj but only to demonstrate the approximately quadratic relationship 
between Ka and the asymmetry parameter. In addition, it should be un- 
derstood that this expression provides only an "average" value for K T and 
the available data is not sensitive to higher-order effects like the "surface" 
part of this term. 

From such an analysis of the ISGMR data in the Sn isotopes, we had 
obtained a value of K T = -550 ± 100 MeV (see Fig. 4 in Ref. [5j). An 
identical analysis of the ISGMR data in the Cd isotopes gives a prelim- 
inary value of K T = —480 ± 100 MeV (see Fig. 2). These numbers 
are in good agreement with each other, and are also consistent with the 
value of K T = — 370 ± 120 MeV obtained from an analysis of the iso- 
topic transport ratios in medium-energy heavy- ion reactions [TU], the value 
K T = — 5001^qq MeV obtained by Centelles et al. [2D] from constraints 
put by neutron-skin data from anti-protonic atoms across the mass table, 
and K T = —500 ± 50 MeV obtained by Sagawa et al. by comparing our 
Sn ISGMR data with calculations using different Skyrme Hamiltonians and 
RMF Lagrangians [21] , A more precise determination of K T will likely re- 
sult from extending the ISGMR measurements to longer isotopic chains. 
This provides strong motivation for measuring the ISGMR strength in un- 
stable nuclei, a focus of current investigations at the new rare isotope beam 
facilities at RIKEN, GANIL, GSI, and NSCL. 

From the ISGMR and ISGDR data, one now has a consistent value for 
Koo = 240 ± 20 MeV [22l [231 El ES] Combined with the "experimental" 
value of K T , this can provide a means of selecting the most appropriate of 
the interactions used in EOS calculations. For example, this combination 
of "experimental" values for K^ and K T appears to rule out a vast ma- 
jority of the Skyrme-type interactions currently in use in nuclear structure 



1 We are using a more conservative uncertainty in the value than that cited previously. 
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Fig. 2. Systematics of the difference Ka — KcouiZ 2 A~ 4 ^ 3 in the Cd isotopes as a 
function of the "asymmetry-parameter" ((N-Z)/A); K Co ui = -5.2±0.7 MeV [21] . 
The solid line represents a least-square quadratic fit to the data. The dashed line 
shows, for comparison, results with the value K T — -550 MeV previously obtained 
from the Sn data. These results should be considered preliminary. 



calculations [261 12Z] • 

Pearson et al. [28] have recently suggested that our data for the Sn 
isotopes may be fitted with K T = -300 MeV, a value they claim is more 
amenable to many theoretical calculations. They fix the K T term at this 
value and by introducing a higher-order term, KsgA 1 ^ 3 , with K T in the 
expansion, obtain a reasonable fit to the data. In effect, they have changed 
the fitting expression from y = A + Bx 2 (which we had used) to y = A' 
+ (B'+C')x 2 and fixed A' and B', to obtain a value for C from fit to 
the data. Clearly, there are several problems with their approach. First, 
considering that C would change very little (~3%) over the range of the fit 
while the x 2 term changes by 230%, C is, for all practical purposes (and 
definitely for the uncertainties involved in the extraction of the coefficient), 
a constant. So, the fit of Pearson et al. replaces one constant with two and, 
of course, one can find a large number of combinations of B' and C that 
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would give essentially the same fit as the one given by B! [Because Pearson 
et al. choose to fix B' at -300 MeV, they obtain a value of -1900A 1 / 3 MeV 
for C. Indeed, if they had chosen the value of -100 MeV for B', they would 
get the same fit for C = -2880A 1 / 3 MeV!]. Second, they do not discuss at 
all the "reasonableness" of the value for K ss (-1900 MeV) that is obtained 
in their fit. Nominally, one expects the higher order terms to be smaller 
than the leading terms. In this case, the higher-order term is nearly 30% 
larger! In this context, we note that the value for K ss that they obtain with 
K T = -550 MeV, leads to the higher-order term being only about 30% of 
the leading-order term, a much more reasonable situation. Finally, there is 
also the issue of their use of a value of 370 MeV for the surface term K s f 
{=cK vo i in Eq. 2 above). Per almost all theoretical estimates, the value of 
c ~ 1 (see Refs. [13 [21]). Pearson et al. obtain c = 1.54 but, again, do not 
justify a number at such variance from other theoretical work. 

4. Summary 

We have measured the strength distributions of the isoscalar giant monopole 
resonance (ISGMR) in the even-A Sn and Cd isotopes via inelastic scatter- 
ing of 386- MeV a particles at extremely forward angles, including 0°. The 
asymmetry-term, K T , in the expression for the nuclear incompressibility has 
been determined to be ~ -500 MeV from the ISGMR data and is found to be 
consistent with a number of indirectly extracted values for this parameter. 
Measurements with the new rare isotope beam facilities would go a long 
way in extending the isotopic chains, thus greatly reducing the uncertainty 
in the extracted value of K T . 
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